Proposed Experiment in Two-Qubit Linear Optical 
Photonic Gates for Maximal Success Rates 



A. Matthew Smith 12 , D. B. Uskov 1 ' 3 , M. Fanto 2 , L. H. Ying 1 , 
and L. Kaplan 1 

department of Physics and Engineering Physics, Tulane University, New Orleans, 
Louisiana 70118,USA 

2 Air Force Research Lab, Information Directorate, 525 Brooks Road, Rome, New 
York 13440, USA 

3 Department of Mathematics and Natural Science, Brescia University, Owensboro, 
Kentucky 42301, USA 

E-mail: amos . matthew . smith . ctrOrl . af . mil 

PACS numbers: 42.50.Ex, 03.67.Lx, 42.50.Dv 

Abstract. Here we propose an experiment in Linear Optical Quantum Computing 
(LOQC) using the framework first developed by Knill, Laflamme, and Milburn. This 
experiment will test the ideas of the authors' previous work on imperfect LOQC gates 
using number-resolving photon detectors. We suggest a relatively simple physical 
apparatus capable of producing CZ gates with controllable fidelity less than 1 and 
success rates higher than the current theoretical maximum (S=2/27) for perfect fidelity. 
These experimental setups are within the reach of many experimental groups and would 
provide an interesting experiment in photonic quantum computing. 
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1. Introduction 

Optics remains one of the most promising methods for quantum information processing 
and computing, due to long photon decoherence times, the ease of photon manipulation, 
and the ability to transmit quantum state information over very large distances. Optical 
states have also been proposed as possible buses between matter qubits [HE]. It 
is therefore desirable to be able to manipulate states at the single-photon level via 
gate implementation. Knill, Laflamme, and Milburn (KLM) significantly advanced the 
prospect of single-photon quantum computing in their seminal paper [3], in which they 
overcame the need for nonlinear interactions by using the inherent nonlinearity of photon 
measurements. In this scheme, the computational system is combined with ancillary 
modes, and the gate operation is performed on the enlarged state space. The ancilla 
modes are measured with photon-number-resolving detectors, leaving the computational 
modes undisturbed and in the desired output state provided the measurement is 
successful. In our previous work [U El E], the authors have shown that a combination 
of analytical and numerical techniques may be used to design optimal linear optical 
transformations implementing two- and three-qubit entangling gates. Here we show 
results for non-ideal gates and suggest an experiment to test them. 

The probabilistic nature of quantum measurement implies a trade-off between the 
success rate of the operation (the probability of obtaining the desired measurement 
outcome for the ancillary modes) and the fidelity (the overlap between the actual 
and desired states of the computational system when the ancilla measurement is 
successful). Previously, solutions were obtained that have the maximum possible 
ancilla measurement success probability given the constraint of perfect fidelity for a 
specified transformation [HE]. In practical implementations, however, the goal of 
perfect fidelity may not always be desirable or even attainable. We have therefore 
generalized our previous techniques to the case of imperfect fidelity, and investigated 
the above-mentioned trade-off between the fidelity and success of the linear optical 
transformations. It was found that for sufficiently small deviations from perfect fidelity, 
a single optimization parameter determines the relationship between fidelity and optimal 
success rate 

In section [2] of this work we briefly describe the theory behind our proposed 
experiment. In section [3] we describe a simple experimental apparatus capable of 
producing CZ gates with varying fidelity and success rates. In section [4] we summarize 
the experimental requirements and provide concluding remarks. 

2. Theory 

The input state to the experiment |\E^° mp ) x |\j/ ancilla } j s a product of the computational 
state containing M c photons in N c modes, and an ancilla state containing M a photons 
in iV a modes. The N c computational modes are those on which the actual gate is 
intended to act. Assuming dual-rail encoding, each qubit is represented by one and 
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only one photon in two computational modes, so we have M c = N c /2. The ancilla state 
may in general be separable, entangled, or an ebit state carrying spatially distributed 
entanglement [7j, though here we propose using only a product state of single-photon 
and zero-photon ancillas, which are relatively simple to produce in an experimental 
setting. 

The linear optical device transforms the creation operator af n ^ associated with 
each input mode i to a sum of creation operators J2j ^ij a j • Here U, which contains 
all physical properties of the device, is an N x N matrix, where N = N c + N a is the 
total number of modes. The total input state may be written as a superposition of 
Fock states \^m) = |^i, n 2 , ■ ■ ■ , n N ), where is the occupation number of the i-th input 
mode, and Y^ n i = M c + M a = M is the total number of photons. The input state is 
transformed as 



N i {N 



* out > = fti* ta > = n niE M out)t ] i°) • 



-1 V n 



We note that Q is a multivariate polynomial of degree M in the elements Uij. 

Once the transformation is complete, a measurement is applied to the N a ancillary 
modes. In the case of a number-resolving photon-counting measurement, (^measured I = 
(fcjv c +i, fcjv c +2, . . . , fcjvl, where ki is the number of photons measured in the z-th mode of 
the ancilla. The resulting transformation of the computational state is a contraction 
quantum map |CT P > = A|^ mP >/Pl^in° mp > II [8J, where A = A(U) is defined by 

i|^ n omp ) = (k Na +x, k Nc+ 2, k N \Cl\*- m ) . (2) 

The linear operator A, which maps computational input states to computational output 
states, contains all the information of relevance to the transformation. 

We define the fidelity as the probability that the desired target gate A TaT has been 
faithfully implemented on the computational modes given a successful measurement of 
the ancilla modes: 

, U |TritiTar,2 

F(A) = l - ^4- , 3 

since TrA Ta,r ^ A Ta,T = 2 Mc for a properly normalized target gate. As we are interested in 
deviations from perfect fidelity, we define 5 = 1 — F as our main parameter [6] . 

We define the success rate of the ancilla measurement to be given by an average 
over all computational input states, 

S(A) =Tr(i t i)/2 Mc ||C/|| 2M (4) 

for general complex U. Note that U need not be unitary, as any matrix can be made 
unitary via the unitary dilation technique by adding vacuum modes [HJ [S]. We also 
note that the Hilbert-Schmidt norm = Tr(A J< A)/2 Mc , used in our definition 

of S, is bounded above by the square of the operator norm, ||A|| 2 = (||A|| Max ) 2 = 

Max((^ n omp |^ t ^l*i„° mp )) ! and below by (||i|| min ) 2 = Min((K mp |^IK mp >), where 
the maximum and minimum are taken over the set of properly normalized input states. 
In the limit F — > 1, II A|| Mm /||yi|| Max — > 1, and all definitions of the success rate coincide. 
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The optimization method we have developed maximizes the success probability S 
for a given target transformation A Tar , for given ancilla resources, and for a given fidelity 
level F < 1. This is mathematically equivalent to unconstrained maximization of the 
function S + F/e in the space of all matrices U, where 1/e is a Lagrange multiplier. 
Here e — > + corresponds to maximizing the success probability while requiring perfect 
fidelity (F = 1). As e is increased, the maximum of S + F/e yields linear optics 
transformations that maximize the success S as a function of the fidelity F. Given one 
transformation U that (locally or globally) maximizes success S for a given fidelity F, e 
may be continuously varied to obtain a one-parameter family of optimal transformations, 
tracing out a curve in success-fidelity space. Note that in general the members of these 
families need not be all unitary, however for some gates of interest, including the CZ 
gate, all members of the family are unitary. 

Figure [T] shows optimal results for the CZ gate. Here each point corresponds to a 
unique unitary mode transformation U. As previously reported we find an interesting 
feature of these unitary matrices. The optimal solution with fidelity F — 1 was found by 
Knill to have a surprising form [9], which we have dubbed the "Knill Form" [5], where 
one mode of each qubit is non- interacting, e.g., in the CZ case U acts as the identity on 
modes 1 and 3 (or equivalently 1&4, 2&3, or 2&4). This form has been found to hold 
for the CZ gate and for the TS Toffoli Sign gate (CNOT and Toffoli respectively are 
equivalent to these up to local rotations). 



CS 0.1 



o.os 




Figure 1. Improved success rates for compromised 5. Numerical results show 
the maximal success rate S vs. S for the CZ gate, and the curve is a best fit to 
the theoretical form Each data point rorrresponds to a distinct linear optical 

transformation, given by a unitary matrix. 



We fit the data for the CZ gate to the known general analytic form [6], 

S(F) = S + S 1 {1- F) 1 ' 2 + S 2 (l — F) + . .. , (5) 
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and truncate to 

S{5) = 0.074 + 0.076 5 1/2 . (6) 

When F — 1, S — and the the success S reduces to the 2/27 ~ 0.074 result found 
by Knill [9], which has been confirmed numerically by the authors [4]. The ratio 
Si/ So = 1.03 ± .03 contains the most interesting information about the system, as 
it is a measure of the relative rate of increase in success as the fidelity is compromised. 
We have used this ratio previously to compare the behavior of several gates, namely 
CNOT, CS(90), and the B gate [5]. All of these could be tested by the experimental 
apparatus proposed here, with some modifications. 



3. Experimental Design 



We now propose an experiment that will test the results shown in figure [T] Reck et al. 
have shown that any discrete N x N unitary transformation U can be implemented as 
a multi-port device consisting only of variable transmittance beamsplitters and phase 
shifters [10J. Their method is a decomposition in which each unitary matrix element 
below the diagonal is transformed into zero by a 2 x 2 rotation matrix embedded in an 
N x N matrix which is otherwise equal to the identity. For example, the 2x2 rotation 

which eliminates the element Un,n-i, takes the form 

. ... \ 



acting on modes N and N — 1 

/ 1 • 



Ti 



N.N—l 



V o 



e r sm(wj 
cos(u) 



(7) 



e l<t> cos(o;) 
— sin(w) / 

The method is recursive and requires one iteration for each pair of modes. Finally, we 
obtain 



U {N)T NjN ^T N , N ^ . . . T 2> iD = I , (8) 

where D is a diagonal matrix of phases. The desired transformation U is then 
decomposable as 

U(N) = D- l T 2 ^l . . . Tjv,at-i • (9) 

Physically, each 2x2 transformation T~j is implemented as a variable transmittance 
beamsplitter with a phase plate on one input mode, while D~ x corresponds physically 
to a phase shift on each output mode [10J. 

Thus a generic two-qubit operation, which needs at least N = 7 modes (N c = 4 
computational modes and N a = 3 ancillas) requires a minimum of 21 beamsplitters and 
28 phase shifters. A controlled unitary gate (N = N c + N a = 4 + 2 = 6) requires at least 
15 beamsplitters and 21 phase shifters. If unitary dilation is required (as is often the 
case) the number of optical elements increases rapidly. However our experiment does 
not require unitary dilation and furthermore as noted by Reck et al., if an element of 
the unitary matrix is already zero, then no transformation is required. The element is 
skipped. 
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Here we return to the "Knill Form," where in the case of CZ we find that nine 
of the elements below the diagonal are already zero. Therefore the unitary transform 
can be implemented with only six beamsplitters and 10 phase shifters. We can perform 
this decomposition for each data point in figure [TJ and find the rotation angles 
and phases 0jj in each case. Surprisingly we find numerically that all of the phase 
shifts, (pi j, are constant along the entire length of the curve in figure [TJ Therefore 
only the six beamsplitter rotation angles Uij out of a total of 36 possible variables 
need to be modified to vary 5, making the experiment much more physically realizable. 
To be specific, the transformation only requires beamsplitters acting on the following 
mode pairs: = (6,5), (6,4), (6,2), (5,4), (5,2), (4,2). Figure || shows that the six 
beamsplitter rotation angles change smoothly with S. Implementing such rotations and 
constant phase shifters will recreate the unitary matrices from figure [TJ 
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Figure 2. Beamsplitter rotation angles. Numerical results showing the six 
beamsplitter rotation angles u>ij vs. 5 for the CZ gate. 



This system lends itself to being implemented with 2x2 Mach-Zehnder 
interferometers (MZI) in place of standard beamsplitters. The transmittance of the 
MZI is controlled dynamically by adjusting the phase difference, without having to 
alter the physical system. These interferometers have already been put on optical 
chips by Thompson et al. [UJ among others. Indeed, significantly larger electro-optical 
matrix switches have been proposed and built for broadband optical communication 
networks [TT1IT2]. 

Figure [3] shows a multi-port device that mixes seven input/output modes (thin lines) 
using 2x2 variable transmittance beamsplitters (rectangles), each of which has a phase 
shifter on one of its input modes (ellipses). An additional phase shifter is placed on each 
device output mode. The thick line is a simple mirror. J. L. O'Brien recently proposed 
a similar 7x7 single-chip MZI-based device made from lithium niobate waveguides [UJ . 
The intended purpose of this chip was to be able to perform any two-qubit unitary 
operation, i.e. any transformation in S77(4). However, such a device would also be 
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~ ~ VariableTransmittance Beamsplitter Optical Mode 

• Phase Shifter Mirror 

Figure 3. General seven-mode optical device for two-qubit gates. A schematic 
depiction of a multi-port optical device capable of performing our experiment as well 
as implementing any generic SU(A) transformation in the measurement-assisted KLM 
scheme [TU] , 

Physical realization of work such as this can occur in many different experimental 
configurations. A bulk optical component configuration allows the greatest flexibility in 
system reconfiguration. The components required to construct this system are readily 
available at fairly low cost. Unfortunately the system suffers from a large and badly 
scaling footprint size, thermal instabilities, and misalignment issues. 

The next logical step would be to implement this system in optical-fiber-based 
components, thus alleviating the need for critical alignment and reducing the overall 
footprint size. Fiber allows for a reconfigurable system with properly connectorized 
components and drastically reduces the setup time for this experimental configuration. 
The fiber-based components required for this experiment are also all readily available. 
The phase shifters can be implemented in a manual stress-induced fiber phase shifter 
or a lithium niobate phase modulator. The variable transmittance beamsplitters can be 
replaced with 2x2 evanescent couplers whose output port transmittance can be changed 
by manual adjustment of the gap between the fibers. The tuning of these devices is slow; 
most are manually tuned and would require simultaneous tuning across the entire circuit. 
Another possible component configuration consists of fiber-coupled bulk beamsplitters 
arranged to form a Mach-Zehnder interferometer with a lithium niobate phase shifter 
inserted into one of the arms. This would allow for high-speed and remote tuning of the 
transmittance. Lastly, the beamsplitter could be implemented with a single device, a 
2x2 lithium niobate MZI. This device allows for the highest-speed reconfiguration and 
transmittance tuning of any of these possible suggested configurations. The respective 
transmittance set for any particular MZI must be monitored due the thermal drift of 
the chip and the effect is amplified due to the modulators sinusoidal transfer function. 
Many commercial MZI's address this problem with thermal stabilization of the chip and 
output power monitoring via an integrated photodetector. 

However, fiber does pose a number of other challenges in a configuration such as 
this. Polarization mode dispersion (PMD) occurring from random imperfections and 
asymmetries can lead to phase variations and these would need to be compensated for 
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with the installation of a polarization controller. PMD can be negated at the cost of 
accepting slightly higher fiber loss with the use of polarization-maintaining (strongly 
birefringent) fibers, in which the polarization is confined to one transmission plane. 
The other related issue to be aware of is polarization-dependent loss, again due to 
asymmetries, where one polarization experiences a higher loss rate. Both of these effects 
can lead to transmittance differences through the interferometers in the experiment. 
Thermal expansion can also induce phase variations in the fiber-based interferometers. 
This can be actively compensated for with the addition of a Pound-Drever-Hall type 
feedback scheme, at the cost of increased system complexity [13]. Lastly, connector 
interfaces between components can lead to photon loss and state degradation. 

The ideal implementation would be a monolithic integrated circuit, which has a 
minimal footprint and provides greatest control of thermal stability. There are a number 
of choices of material (silica, silicon, lithium niobate, etc.), but the final decision depends 
on the wavelength of interest, photon detection capabilities, and the intended function 
of the circuit. For example, silica will have lowest loss for passive waveguides in the 
telecom regime, but active devices in silica must be thermal and are therefore slow [14J. 
Lithium niobate will have higher loss waveguides than silica in the telecom regime, but it 
is far superior for active components such as modulators due to its nonlinear properties. 

Initial circuits including photon sources, waveguides, and detectors may have to be 
built on multiple chips, each with its own respective components that are best suited for 
being mated together and then eventually merged onto one monolithic chip. Regardless 
of the chosen material, integrated waveguides allow movement towards scalable photon 
circuits [TT] . 

Photon detection is the final piece of the puzzle. Standard avalanche photodiodes 
(APDs) used with the detectors at 1550 nm have a higher noise than those at 800 
nm, and neither are capable of photon-number resolution. A second option is to use 
superconducting single-photon detectors, e.g. niobium nitride nanowire detectors [15J, 
which can be configured to be capable of photon number resolution. The third option 
is to employ transition-edge detectors, which have demonstrated some of the highest 
number-resolving system detection efficiencies to date, at 98%. 

4. Conclusion 

We have shown the theoretical basis and interest for this experiment. At this time it 
is the only apparent means of experimentally confirming the key ratio Si/Sq, which 
quantifies the trade-off between fidelity and success, for the CZ or CNOT gate. The 
experimental setup may naturally be extended to explore the behavior of other quantum 
gates of interest. The components needed for the execution of the experiment are well 
within the means of many experimental groups. The main stumbling block is the expense 
of purchasing number-resolving detectors. However, any group already possessing these 
detectors should be able to implement this scheme with relative ease. 



Proposed Experiment in Two-Qubit Linear Optical Photonic Gates 



9 



Acknowledgments 

This work was supported in part by the NSF under Grants PHY-0551164 and PHY- 
0545390, and by the National Research Council Research Associateship program at the 
AFRL Rome Research Site, Information Directorate. 

References 

[1] Moehring D L et al. 2007 Nature 449 41 
[2] Kimble H J 2008 Nature 453 1023 

[3] Knill E, Laflamme R, and Milburn G J 2001 Nature 409 46 

[4] Uskov D B, Smith A M, and Kaplan L 2010 Phys. Rev. A 81 012303 

[5] Uskov D B, Kaplan L, Smith A M, Huver S D, and Dowling J P 2009 Phys. Rev. A 79 042326 
[6] Smith A M, Uskov D B, Ying L H, and Kaplan L 2011 Phys. Rev. A 84 032341 
[7] Wilde M M and Uskov D B 2009 Phys. Rev. A 79 022305 

[8] Kraus K 1983 Lecture Notes: States, Effects and Operations: Fundamental Notions of Quantum 
Theory (Springer, New York) 

[9] Knill E 2002 Phys. Rev. A 66 052306 
[10] Reck M, Zeilinger A, Bernstein H J, and Bcrtani P 1994 Phys. Rev. Lett. 73 58 
[11] Thompson M G, Politi A, Matthews J C F, and O'Brien J L 2011 LET Circuits Devices Syst. 5 94 
[12] Sohma S et al. 2006 Euro. Conference on Optical Communication pp 1-2 
[13] Drcver RWd al. 1983 Appl. Phys. B 31 97 

[14] Matthews J C F, Politi A, Stcfanov A B, and O'Brien J L 2009 Nature Photonics 3 346 
[15] Natarajan C M et al. 2010 Appl. Phys. Lett. 96 211101-1 

[16] Nielsen M A and Chuang I L 2000 Quantum Computation and Quantum Information (Cambridge 

University Press, Cambridge) 
[17] Lapaire G G, Kok P, Dowling J P, and Sipe J E 2003 Phys. Rev. A 68 042314 



